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Abstract. The presence of Majorana zero-energy modes at vortex cores in a 
topological superconductor implies that each vortex carries an extra entropy so, 
Cu i given by (fca/2) In 2, that is independent of temperature. By utilizing this special 

property of Majorana modes, the edges of a topological superconductor can be 

cooled (or heated) by the motion of the vortices across the edges. As vortices 

T~j ' flow in the transverse direction with respect to an external imposed supercurrent, 

(-H | due to the Lorcntz force, a thermoelectric effect analogous to the Ettingshausen 

^\ , effect is expected to occur between opposing edges. We propose an experiment to 

f\ ■ observe this thermoelectric effect, which could directly probe the intrinsic entropy 

of Majorana zero-energy modes. 
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£\j i 1. Introduction 



The search for Majorana modes in condensed matter — a subject of intense 
experimental effort — is driven in large part by the expectation that whenever such 
fcrmions appear as zero-energy modes bound to either vortices [H [2] or end points 
of superconducting quantum wires [3J 21 [5] they are characterized by non-Abelian 
braiding statistics [6j[7l[8]. Such particles could, as a result, be utilized for quantum 
information processing [5J. 

A number of condensed matter systems, most notably the v = 5/2 fractional 
quantum Hall state [TU] and chiral p-wave superconductors [TJ [2] , are expected to host 
such non-Abelian quasiparticles. Alternatively, chiral superconductors supporting 
Majorana zero-energy modes can be fabricated as heterostructures composed of an s- 
wave superconductor and either a topological insulator [11] or a semiconductor having 
strong spin-orbit coupling and an additional source of Zeeman splitting [12l [T3l E] ■ A 
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modification of the latter scheme for the case of a semi-metal may remove the Zeeman 
splitting requirement |15j . 

However, an unambiguous experimental observation of a Majorana zero-energy 
mode remains elusive, thus far. Following earlier proposals for the interfcrometric 
detection of non-Abelian anyons in the v = 5/2 fractional quantum Hall state [16j [T7l 
I18j , similar ideas were put forward in the context of topological superconductivity |19| 
EOlET]. Another possible signature of Majorana modes would manifest itself through 
an unusual 47r (rather than the conventional 2tt) periodicity of a Josephson current as 
a function of the phase difference across the junction [3J[52]. A zero-bias tunnelling 
anomaly and the corresponding 2e 2 jh quantization of the tunnelling conductance from 
a single metallic channel into a Majorana zero-energy mode (23] [24] could provide 
another signature. While the first experimental results consistent with the latter 
prediction have been just reported [25j . they cannot not address the most interesting 
feature of these quasiparticlcs, namely their non-Abelian statistics. 

This brings us to another unique but less explored feature - an intrinsic zero- 
temperature entropy of so = (fcB/2)ln2 per Majorana zero mode - which is 
also an essential hallmark of such quasiparticlcs. This entropy results from the 
exponential growth of the ground-state degeneracy with the number of quasiparticles, a 
precondition for their non-Abelian statistics [9]. Hence, a measurement of the intrinsic 
entropy carried by each vortex can be taken not only as an unmistakable signature 
of Majorana zero modes but also as an indication of their unusual statistics. It has 
been argued by Yang and Halperin that the presence of this zero-temperature entropy 
leads to an enhancement of thermopower j26| . Furthermore, it can be utilized for the 
adiabatic cooling of systems supporting non-Abelian anyons (27] [28] . 

In this paper, we show that the zero-temperature entropy carried by vortices 
in a topological superconductor induces a magneto-thermoelectric effect - the 
Ettingshausen effect. We also propose a specific setup — a heterostructure combining a 
topological insulator and an s-wave superconductor with a wide Josephson junction — 
that should give rise to a measurable signal for this effect under plausibly realizable 
experimental conditions. 

2. Thermoelectric effect 

Let us begin by introducing an intuitive qualitative picture of the edge thermoelectric 
effect in a 2D topological superconductor with broken time reversal symmetry. An 
edge of such a superconductor is characterized by the existence of a gapless chiral 
mode. (Depending on the net vorticity inside the topological superconducting region, 
such a mode may or may not be exactly at zero energy.) When a vortex enters the 
topological region, it necessarily crosses this gapless edge and changes its spectrum; 
as a result a quantum state is "peeled off" from the edge to form the Majorana zero- 
energy mode localized at the vortex core. This process, in turn, reduces the entropy 
associated with the edge modes by exactly so = ^&Bln2, which is carried away by 
the vortex. In the reverse process, whereby a vortex moves out from the topological 
superconductor region, the same amount of entropy is added back into the opposing 
edge. 

Alternatively, we could analyze the effect of a pair of vortices entering into the 
topological superconductor. The advantage of this approach is that we need not 
concern ourselves with the details of the edge spectrum reconstruction upon the 
passage of each successive vortex. Once the pair has moved deep into the topological 
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region, the energy spectrum of the edge must return to its original form. On the 
other hand, both vortices now carry zero modes (provided they are sufficiently well 
separated). Assuming that the temperature is lower than the dimensional quantization 
energy inside their cores (the minigap) , there are no other entropy changes associated 
with the vortices. On the other hand, if the vortex passage is an adiabatic process, the 
total entropy associated with the vortex and edge states must be conserved. The only 
way to reduce the entropy of the edge states in order to compensate for the entropy 
now carried by the two vortex zero modes is by reducing the temperature of the edge. 
(Only the those states that are within the 'temperature window' of the edge contribute 
to its entropy.) Naturally, the opposite effect results from the pair of vortices leaving 
the topological region. 

Owing to the Lorentz force, vortices may be driven across the sample via an 
externally imposed supcrcurrent. Because such vortex flow cools one and heats 
the other of an opposing pairs of edges, a thermoelectric effect, analogous to the 
Ettingshausen effect, occurs between opposing edges. As we shown below, such 
an effect can be quantified in terms of the ratio of the temperature difference AT 
between the opposing edges and the voltage drop V in the applied current direction 
AT/V = (e/fc B )(12/7r 2 )ln2. 

As discussed by Yang and Halperin [26] , the presence of the intrinsic entropy sq 
per vortex is justified only when the temperature is higher than the energy splitting 
of the zero-energy modes, i.e., T ;§> To ~ Ae~ l ' l °. For the present work, A is the 
superconductor gap, lo is the typical size of vortex, and / is the distance between 
vortices. Thus, in the limit of dilute vortices, To is exponentially suppressed. At a 
nonzero temperature these vortices carry additional entropy, due to other minigap 
states |29| . However, if the temperature is lower than the minigap, i.e., k^T <C £ m g, 
these contributions are suppressed exponentially by S mg ~ |£ mg /T|e~l' £ms /' SBT l and 
can be simply ignored (see | Appendix A[ ). Therefore, the edge thermoelectric effect 
due to the presence of the intrinsic entropy will only be prominent in the temperature 
range T <^T <C £ mg /k B . 

As a concrete example, we shall consider the schematic setup shown in Fig. [TJ 
which illustrates the mechanism underlying this thermoelectric effect. We consider 
a topological insulator (TI) that interfaces with an s-wave superconductor, so that 
a superconducting pair potential is induced in the contact region of the TI via the 
proximity effect. This region effectively emulates a 2D topological superconductor. 
Each vortex in this region would have a Majorana zero-energy mode bound at its 
core [IT] . We envision magnetically gapping the rest of the TI surface by depositing 
on it a ferromagnetic insulator. (We shall revisit this point and discuss more practical 
means of producing such a gap later on.) In this setup, a one-dimensional chiral 
Majorana edge state will form at the boundary of the superconducting region P3B [5D] . 
We then imagine subjecting this region to a transverse magnetic field of strength 
B > B c i, which will result in a vortex density n v = -B/$o = B/(h/2e) in the 
superconducting slab. Here, B cl is the first critical field and $o = h/2e is a 
superconducting flux quantum. Finally, we envision applying an external current 
to the superconductor, which will induce vortices to move laterally, between the two 
opposing edges. 

If the vortices move with velocity u (which depends on the frictional force on 
moving vortices) , the entropy current in the TI would be given by |30[ 131] 

2eB 
is = s n v u = s —— u. (1) 

h 
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Figure 1. Schematic plot of the proposed setup showing the thermoelectric 
mechanism for the edge state. A 3D topological insulator (colored yellow) is in 
contact with an s-wave superconductor (colored grey) from one side, and the rest 
of surface is coated by ferromagnetic insulator of magnetization M normal to 
and inward from the surface indicated by the green arrow at each surface. A 
circulating Majorana edge state (shown in red) forms at the boundary separating 
the regions of superconducting pairing and magnetic gaps. A magnetic field is 
applied in the 2-dircction. A current driven in the y-direction produces a force 
that pushes vortices in the x-direction. The passage of a vortex cools the left edge 
and heats the edge on the opposite (i.e., right) side. 



In order to sustain a constant vortex motion, a uniform electric field E = B x u should 
be applied in the direction perpendicular to both the magnetic field and the vortex 
motion. Hence, in terms of the applied electric field, the entropy current becomes 



2e 
Js = s o-^" E x B, 



(2) 



where B = B/ |B| is the unit vector in the direction of the magnetic field. Here, the 
effect of the Magnus force has been ignored, as it will not affect the conclusion given 
in Eq. ©. 

To understand the heating (cooling) of edges due to the vortex flow, we first need 
the heat capacity per unit length C v of the chiral Majorana edge state, which is given 

by 

c v = -jklTp EFl (3) 

where pe f = l/(4-7r?iv^) is the density of states at the Fermi energy, and v^ is the 
velocity of the edge states. The energy current flowing out of/into the heated/cooled 
region is given approximately by 

— — = c v v^,dl = 

dt v v I2h 

where 5T is the temperature variation due to heating/cooling. By balancing this 
energy flow with the heat added to (or removed from) the edge states due to the 
vortex motion crossing the edge, i.e., 



-klTST, 



(4) 



dQ v 
~dT 



LTjs = Ts 



2eV 



(5) 
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we obtain the result 

ST = — —T~ ' ( 6 ) 

where L is the length of the heated region and V = L\E\ is the voltage drop across 
the superconductor. 

As the amount of entropy removed from one edge is deposited by vortices at the 
other edge, the ratio of the temperature difference AT between the opposing edges 
and the voltage drop is given by 

AT _ W _ 12M2 . m 10< 

V V it* Kb 

Because both the temperature difference and the voltage drop directly originate from 
the motion of vortices, any vortex pinning should not affect this signal [32] • As this 
thermoelectric response is quite substantial (i.e., AT rs 10 mK for 1 pN of applied 
voltage), it should be possible to measure this effect provided it proves possible 
to measure the edge state temperature while keeping the edges isolated from the 
environment. 

Although the setup in Fig. [1] is useful for demonstrating the idea of the 
thermoelectric effect at a conceptual level, we emphasize that in reality the effect 
can be masked, and therefore difficult to measure, in this simple setting, due to the 
following reasons. Firstly, an Abrikosov vortex in an s-wave superconductor possesses 
a normal core, and thus carries entropy in addition to the contributions from the zero 
mode and minigap states. We note that these additional sources of entropy, although 
oblivious to the existence of the edge states, can build up a temperature gradient — -the 
classical Ettingshausen effect — in the bulk of the superconductor, and hence obscure 
the edge thermoelectric effect [3T] [32] . 

Secondly, the motion of vortices may not be strictly perpendicular to the direction 
of the applied current, instead having a Hall angle induced by the Magnus force and 
depending on materials properties [3T] [33] . 

Thirdly, the motion of a vortex in a superconductor induces a non-zero resistivity 
that is proportional to the applied magnetic field, p oc p n B/B C 2, where p n is the 
normal-state resistivity and B c2 is the second critical field strength [30]. Because 
conventional superconductor materials are characterized by a small normal-state 
resistivity, a superconductor with moving vortices yields a small resistivity. For a 
fixed voltage, the smaller the resistivity the larger the Joule heating. As a result, the 
superconductor can be heated considerably due to the motion of vortices. 

Finally, carriers in chiral edge states in a spatially extended heating (or cooling) 
region may equilibrate with the environment before they can reach the thermometers. 
This could make it difficult to measure the temperature difference between two edge 
states that results from the transfer of the intrinsic entropy carried by the vortices. 

3. Wide Josephson Junction Device 

In order to overcome the aforementioned obstacles to detecting the edge thermoelectric 
effect, we now propose an alternative device that utilizes Josephson junctions, as shown 
in Fig. [2] A wide Josephson junction, in which Josephson vortices can propagate, is 
situated under a slab of TI. A constant supercurrent density J s is applied across 
the junction (i.e., in the y-direction), in order to push Josephson vortices along the 
junction (i.e., in the x-direction) by means of a Lorentz force. An impedance-matched 
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Figure 2. Schematic of a setup involving a wide Joscphson junction. 
An insulating tunnel barrier (white region) sandwiched between two s-wave 
superconductors (gray regions) forms a wide Josephson junction. The geometry 
of the junction is as follows: the thickness of insulator is I; the depth of junction 
area is d; the width is w. The width of the TI is denoted by w TI . The lighter 
and darker blue ellipses, located where edge states cross the Josephson junction, 
indicate cooling and heating regions, respectively. Josephson vortices arc formed 
by applying a weak magnetic field B z in the z-direction. The magnetic fields 
are also applied, indicated by green arrows, in order to gap the surface states 
in both positive x- and negative y-directions. The chiral Majorana edge states 
appear between the superconducting and magnetic gapped regions as indicated 
by red perimeter with the red arrows indicating the flow directions of Majorana 
edge states. Between magnetic gapped surface states with opposite magnetization 
directions, additional edge states appear and connected to the Majorana edge 
states as indicated by two vertical red arrows. Here, the upper surface state of the 
topological insulator is not gapped as B z is very weak and is in general screened 
by the superconductor. An impedance-matched resistance circuit is attached at 
one end of junction in order to guide the vortex out of the junction. Finally, a 
supcrcurrent J B is applied through the junction to drive the vortices. 



resistance circuit should be placed at one end of the junction, so that vortices are not 
reflect from the edge of the junction [MJ |3S]. Finally, to induce a magnetic gap on 
the free surface of the TI, one could apply a magnetic field parallel to the interface; 
a chiral Majorana edge states then form at the boundary of the superconductor, as 
shown in Fig. [2] 

As with an Abrikosov vortex, a Josephson vortex carries a Majorana zero mode 
in the topological superconductor region [21] . Therefore, an analogous thermoelectric 
effect, caused by the mechanism described in the previous section, should occurs 
when dilute Josephson vortices move across the sample. A dissipativc motion of 
such a vortex along the junction induces a voltage pulse across (due to the Josephson 
relation) ; the time-averaged voltage drop across the junction is then given by 

V = ^$o, (8) 

where v is the passage frequency of Josephson vortices through the junction. One 
can subsequently show that the chiral edges emerge downstream from the junction, 
heated or cooled according to Eq. ((6j) with a simple substitution V h-> V. As 
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the cooling/heating processes only take place at or near the junction, it should be 
possible to probe the thermoelectric effect before the edge state equilibrates with the 
environment. 

The advantage of dealing with wide Josephson junctions can be understood 
qualitatively, before providing a formal treatment of its soliton excitations (i.e., 
Josephson vortices). Although the superconducting phase winds by 2tt around both 
types of vortices, a Josephson vortex has a phase core but no minigap states, whereas 
an Abrikosov vortex contains a normal core with minigap states [36]. Thus, the 
propagation of Josephson vortices in a conventional superconductor carries no entropy, 
resulting in a small temperature gradient, at most. Furthermore, as the friction and 
pinning forces encountered by a Josephson vortex in a well-fabricated junction can be 
much smaller than those associated with an Abrikosov vortex moving inside a bulk 
superconductor, we expect much less heat dissipation due to the motion of a Josephson 
vortex. 

The dynamics of a wide Josephson junction can be described by a sine-Gordon 
equation that includes damping and driving forces |37j , cf . |Appendix B| 

d 2 d 2 d\ 

0^2 ~d^2 - a fr) ^(C ) r)=sin^(C,T)+ 7) (9) 

where <p(£, r) is the position and time-dependent gauge-invariant phase difference 
across the Josephson junction as a function of dimcnsionless variables £ = x/Xj and 
t = ct/Xj. Here, 



2irno(2X L +£)J c : ~ V e Mo(2A L + I) ' 



*J = W n- /ox ■ ,W . 5 =\l - ,m I « > ( 10 ) 



are the Josephson penetration depth and the effective speed of light, respectively. 
These parameters have the meaning of a characteristic size scale and a propagation 
speed of a Josephson vortex along the junction; they are governed by material 
properties and geometric parameters of the Josephson junction: the permittivity e of 
the insulator, its thickness £, the London penetration depth of the superconductor Al, 
and the critical current density J c of the Josephson junction. The damping coefficient 
a = /.jo(2Al + £)cXj/£p is inversely proportional to the resistivity p of the insulator, 
and the (dimcnsionless) driving force 7 = J s / J c is the ratio of the supercurrent and 
critical current densities. 

To solve the sine-Gordon equation, one also needs to specify the boundary 
conditions. For a junction of width £ , tli e boundary condition for Eq. ^ reads 

^(Co)-^(0) = ^B 2 (2A l +^)AjCo, (11) 

in the presence of applied magnetic field B z . Indeed, this condition simply states 
that the total phase winding along the junction has to match the total flux threading 
through the junction. 

In the limit a, 7 <C 1, one can first ignore the contributions from the damping and 
driving forces. When the total magnetic flux in the junction area is exactly one flux 
quantum, the sine-Gordon equation (j9|) has a soliton solution having a profile given 
byEB 



tp(x, f) = 4 tan 



/ x — vt 

exp ±- 



'x,^i-(v/cy 



(12) 
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in terms of the real time and coordinate. This moving soliton, trapping exactly one 
flux quantum within the size If = Ajyl — (v/c) 2 , is a Josephson vortex. Here, the 
propagation velocity v of the Josephson vortex is determined by the balance between 
the damping and driving forces and, takes the value |37j 

v = c/y / l + (4a/ 7 r 7 ) 2 , (13) 

which can be controlled by the ratio of damping coefficient and driving constant. In 
the dilute limit (i.e., vortex density d v < 1/Aj), the phase profile at the Josephson 
junction increases monotonically, and is roughly a train of isolated Josephson vortices 
moving with velocity v. As a result, the vortex density d v matches the magnetic 
flux density threading through the junction, and hence can be adjusted simply by 
controlling the magnetic field strength B z . 

From energy conservation, we have that the heat dissipated per unit length by 
a propagating vortex is precisely the work done by the Lorentz force acting on the 
vortex, i.e., 

P = $o J s v. (14) 

To estimate how much energy is transferred to the edge state due to this heat 
dissipation, we first assume that the cross-section of the edge state is of order £ s £, 
where £ s and £ are the coherence length of the bulk s-wave and the topological 
superconductors, respectively. Here, £ s provides the penetration depth of the edge 
state in the z-direction, and £ is roughly the size of the edge state in the x-direction. 
Because it takes a time £/v for a vortex to pass the edge state, we estimate that the 
edge state will be heated with energy 

Q = U^oJs- (15) 

which is independent of the propagating velocity. Similarly the total heat dissipation 
for transferring a vortex through a junction of depth d and width w can be estimated 
as Qt = dw&oJ s . Indeed, the heat dissipated by a propagating Abrikosov vortex 
driven by an applied supercurrent is exactly the same as that of a Josephson vortex, 
as given in Eq. (fl4|). However, because it requires a much smaller supercurrent to 
drive the vortex moving along the Josephson junction, heat dissipation becomes a 
much less severe issue for the wide Josephson junction device. 

4. Possible Experimental Realization 

To achieve an appreciable temperature difference, the average voltage V should be 
at least in the range of 0.1 ~ 1 yuV, which corresponds to a passage frequency 
of v = 50 ~ 500 MHz. This voltage would result in a temperature difference of 
AT ps 1 ~ 10 mK between the two opposing edges. With a fixed vortex velocity, 
the passage frequency, and hence V, increases with increasing applied magnetic field 
B z (i.e., the vortex density). Hence, the temperature difference between two edges 
due to the Ettingshausen effect can be tuned by the magnetic field strength. To 
understand the issue of feasibility, let us now show that an Ettingshausen effect having 
a measurably large signal can be established in a wide Josephson junction device within 
reasonable materials parameters. 

In our analysis, we assume the following wide Josephson junction geometry (see 
Fig. [2J: the thickness of the insulator £ = 2 nm, the depth of the junction d = 5 /jm, 
and the width of the Josephson junction w = 0.1 m. As a concrete example, the 
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Josephson junction is constructed by an Al^.0^ insulating layer of e re 10eo sandwiched 
between a pair of s-wave superconductors made of Nb-Sn, and having materials 
properties: the superconductor pairing potential A re 3.4 meV, the coherence length 
£ s = 3.6 nm, and the London penetration depth Al = 124 nm. A Josephson junction of 
this type can be fabricated [34j [35] , and is expected to have a critical current density 
J c ranging from 10 5 to 10 7 A/m 2 . In the present discussion, we assume J c = 10 6 
A/m 2 . 

Using these materials properties, together with the flux quantum <f>o = 2.07 x 
10~ 15 V-sec, from Eq. (fl"Tj| one immediately obtains 

Aj re 32 nm, c re 8.5 x 10 6 m/s. (16) 

With the pairing potential A re 3.4 meV and the critical current density J c = 10 6 
A/m 2 , the tunnelling resistivity of the Josephson junction can be estimated to be 
p re 2 O-m, which corresponds to a damping coefficient a re 0.02 |38j . 

As the vortex propagation speed has to be slower than the Fermi (edge-state) 
velocity of the topological insulator (vf ~5x 10 5 m/s), the wide Josephson junction 
should be operated in the damping-dominated regime 7 <C a < 1. Hence, the size of a 
propagating Josephson vortex is If re Aj since u<c and the width of the topological 
insulator should obey wri 3> Aj. By requiring the vortex velocity vj re 5 x 10 4 m/s, we 
obtain from Eq. (fTB"]) the corresponding supercurrent density of 7 = J s /J c re 1.5 x 10 -4 . 

From Eq. (fTTj) . the density of Josephson vortices is given by d v = B z (2\ L +£)/<& Q . 
As the mechanism leading to the Ettingshausen effect requires a dilute vortex density, 
i.e., d v < 1/Aj, we obtain the constraint B z < 0.25 niT. Then, by using the passage 
frequency v = vd v and the average voltage V = v<&o, we obtain the relation 

V = v(2X h + i)B z . (17) 

In addition, by using the aforementioned parameters, the average voltage is restricted 
to V < 3 /iV, and can be controlled directly by changing the magnetic field B z . From 
Eq. (UJ and Eq. (jTTJ) , we note that temperature difference is AT re 1 ~ 10 niK when 
tuning the magnetic field in the range of B z re 0.01 ~ 0.1 mT. To enhance the strength 
of the Ettingshausen signal, one could introduce a set of wide Josephson junctions, 
situated in parallel and separated by a distance larger than London penetration depth 
Al- The temperature difference AT would then scale linearly with the number of 
junctions. 

A superconductor in contact with a topological insulator not only produces 
superconductivity via proximity but also renormalizes the Fermi velocity of the surface 
state [2H] • H the tunnelling rate between the superconductor and topological insulator 
is optimized, and the chemical potential of the topological insulator is close to the 
Dirac point, the induced pair-potential Ati re A/2 is about half of the bulk s-wave 
superconductor pair-potential, and the Majorana edge-state velocity v^, re vp/2, i.e., 
the renormalized Fermi velocity is about half original Fermi velocity [40]. Thus, the 
coherence length of the induced topological superconductor is £ = hv^/A re 25 nm. 
Then, by using aforementioned material parameters, we find that the minigap of a 
Josephson vortex in the topological superconductor region is estimated to be 



^ g «^^«0.13meV, (18) 

which equals to 1.5 K, cf. |Appcndix C[ Taking for the operating temperature T ~ 0.1 
K re (£^ /fca)/15 K, which should be readily achievable experimentally, the entropy 
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contribution from minigap states is about S mg ~ 0.0014 s an d can be neglected, 
(cf. |Appcndix A[ ) . 

Compared with the setup propagating Abrikosov vortices, the issue of heat 
dissipation is dramatically lessened for the wide Josephson junction setup. From 
Eq. ((i"5j) , we see that the heat transferred to the edge state due to a propagating 
Josephson vortex can be estimated as Q w 10~ 28 J. Thus, compared with the 
heat added to (or moved from) the edge states due to the crossing of a vortex, 
(fceT/2)ln2 ~ 5x 10~ 25 J, we can ignore heating due to dissipation from vortex 
motion. In addition, the total heat dissipation for a vortex moving through the 
junction is about Qt ~ 10 -18 J/vortex, which should be drained away from the system 
in order to keep it at a constant temperature. 

Finally, we mention that by measuring the thermopowcr voltage AV and the 
conductance of point-contact tunnelling into the edge state, the temperature of 
edge states can, in principle, be inferred via the Mott relation [U 22]. For the 
tunnelling that occurs between two normal metals (or charged ID channels), the Mott 
relation reads S = AV/AT = -7r 2 fc|T(dlnGy<9^)/(3e), where AT is the temperature 
difference between the two metals, G is the conductance, and fj, is the chemical 
potential. Although we are concerned with the tunnelling into a charge-neutral edge 
state in a superconductor system, we expect that the Mott relation should hold, up 
to an overall prefactor. We note that establishing a Mott relation for superconductors 
is in itself an interesting question that is worth careful examination in the future 
work. We should also emphasize that probing the temperatures of edge states can 
be a challenging experimental task. However, we envision that a setup akin to the 
measurement of quantum Hall edge states, as in Ref. [42], could be a viable scheme. 

Conclusion 

In summary, we have shown that the intrinsic entropy carried by vortices possessing 
Majorana zero-energy mode leads to an Ettingshausen effect between opposing 
sides of Majorana-carrying edge states. In addition, we proposed that this effect 
could be measured using a wide Josephson junction situated on a superconductor- 
topological insulator heterostructure, and we have shown that this setup should permit 
measurement of the Ettingshausen effect within the range of experimentally accessible 
parameters. At low temperature, i.e., T <C £mg/&B; this unique thermoelectric effect 
can be related to the intrinsic entropy, and thus provides a distinct probe of the 
non-Abelian nature of Majorana fermions. Moreover, this edge Ettingshausen effect 
could potentially be used as a refrigeration process for cooling small objects, such as 
a quantum dot. 
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Appendix A. Entropy contribution from the minigap states 

When a vortex moves from a normal superconductor into a topological superconductor 
region, it acquires a zero-energy quantum state. The entropy per vortex associated 
with such a state is sq = (feB/2)ln2. In principle, the energy levels associated with 
other minigap states may also be affected when a vortex is moving from one region to 
another, in which case there will be an additional contribution to the heat transport. 
The entropy carried by a fermionic state with energy £ j at temperature T is given 
by 

Si = -k B { Vl In(pi) + (1 - Pi) ln(l - Pl )} , (A.l) 

where pi = l/(exp(£j/fcBT) + 1) is the Fermi-Dirac distribution function for chemical 
potential fi = 0. Because the entropy is extensive, the total entropy of multiple 
quantum states can be simply added, to give 

£ng = l>Si, (A.2) 

i 

where gt is the degeneracy. Because the entropy is suppressed exponentially by a 
factor exp(— \£i/k-eT\) in the limit \£i/k&T\ 3> 1, the lowest energy state (with energy 
£mg) of the minigap states makes the most substantial contribution, and the entropy 
can be approximated by 

En.. 

exp 



Sme ~ &R 1 



'nig 



knT 



knT 



(A.3) 



Therefore, one can conclude that the entropy contribution from the minigap states 
can be ignored at sufficient low temperatures. 

Appendix B. Sine-Gordon equation 

Here, we derive the equation of motion for a wide Josephson junction in terms of the 
phase difference across the junction. Our derivation emphasizes the relation between 
the magnetic field distribution at the junction and its resulting phase, magnetic field, 
and current distributions inside the bulk superconductors. 

Referring to the coordinates of Fig. [2] let us consider a Josephson junction 
composed of an insulating layer sandwiched between two semi-infinite superconductors 
at areas y > £/2 and y < —1/2. Hence, any z dependence can be ignored. The 
magnetic field B z (x,t) penetrating through the insulating area, —£/2 < y < £/2, is 
assumed to have no y dependence, and thus obeys the Maxwell equation 

f)Fl / 1 f)E \ 

Ho[j c sm<p + J s + -E y + e—^), (B.l) 



dx f-"y"-~—r ' - ■ p ~v • dt 

where ip(x, t) is the gauge-invariant phase difference between two superconductors, /xo, 
e and p are, respectively, the permeability, permittivity and resistivity of the insulator, 
and E y is the electric field in the y-direction inside the insulator. On the right hand 
side, the first term is the Josephson supercurrent, the second term J s is the bias 
supcrcurrent, the third term is the tunnelling current through the junction and the 
forth term is the displacement current. By using the Josephson relation [31] 

dip _ 2e£ . 
~dt ~ ~fr 



-E y , (B.2) 
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the Maxwell equation becomes 

dB z ( * caV *o d<p\ m _. 

-Ho[J B Bm(p + J B + _ . _ „ + — - - — . (B.3) 



dx ^' v 2tt£ at 2 2^ at 

Now, our goal is to associate B z with tp, and thus derive an equation of motion only 
in terms of ip(x,t). 

For simplicity, we assume that the superconductor coherence length £ vanishes, 
and the healing effect due to the presence of the boundary and the magnetic field can 
be ignored. Because the magnetic field B z at the junction effectively provides the 
boundary conditions for both sectors of superconductor, the magnetic field and phase 
distributions inside the superconductors can be solved as a boundary value problem. 
To further simplify matters, we take the limit £ — > when solving this boundary value 
problem. Also, we only consider the boundary value problem for the y > regions 
and then make use of the mirror symmetry to obtain the solution for the other half- 
plane. 

To derive the magnetic field and phase distributions inside the superconductor, 
we first recall the supercurrent density 

where the prefactor involves the supcrfiuid density n s and the effective mass m*, 
<f) is the superconductor phase, and A is the vector potential. When ignoring the 
capacitance (as we shall do), current conservation, V ■ 3f = 0, leads to 

V 2 0-J^V-A = O. (B.5) 

From V x B = ^3f and B = V x A, we then obtain the Maxwell-current relation 

- V 2 A + V(V • A) = Mo ^ (V0 - ^a) . (B.6) 

To proceed, it is useful to choose a gauge. In order to read out the phase difference 
across the junction, it is most convenient to choose the following one: 

4>{x,V,t) = -<f>(x,-y, £), A x (x,y,t) = -A x (x,-y,t), A y = 0, (B.7) 

in which case the phase difference is immediately given by f(x) = 2(f>(x,y — 1/2). In 
this gauge, current conservation (|B.5P reads 

W-f^ = 0, (B.8) 

and Eq. (|B.6[) explicitly becomes 

d 2 A x 2e/ion s h ( dip 27r 



dy 2 to* \dx $ 



J o 



d 2 A x 2e/i n s hdcj) 

(B-9) 



dxdy to* dy 

Thus by using the relation B z = —d y A Xl we arrive at result 

1 

A 



V 2 B Z - —B z = 0, (B.10) 



V L 



together with the boundary condition B z (y = 0) = B z . Here, Al = \J 'to* / '(2e) 2 [ion s 
is the London penetration depth. 
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30 40 y 50 60 70 



Figure Bl. The solution ip(C) of the sine-Gordon equation is plotted as a function 
of ( for a magnetic flux density <E>o/10Aj. As £ is in units of Aj, vortices are then 
separated by a distance lOAj in (real) space. 



The differential equation (|B.10|) can be solved via the Green function method. 
The Green function with homogeneous Dirichlct boundary conditions for y > is 
given by 

e ik m (x-x') f e ik y {y-y') _ Q ik v (y+y') 



G(x,y;x',y') = - 



d 2 /fc 



(B.ll) 



(2tt) 2 kl + kl + l/\l 

By using Green theorem, we thus have that the field distribution in the upper half- 
plane is given by 

B z (x,y>0) = B(0)c- y/XL + —2\B z {k)\e~ akV cos{kx + k ), (B.12) 

Jo+ 2tt 

where B z (k, f) = \B z {k 1 f)\e l6k is the Fourier amplitude of the magnetic field B z (x, f) at 
the junction, and a k = \/k 2 + l/A^- We have also used the relation B(—k) = B*(k), 
as B(x,t) is real. The solution for y < can be inferred using mirror symmetry. 

Next, observe that the component B z (0) is the averaged magnetic field 
distribution, and that B z {k) captures any spatially non- uniformity. For a uniform 
applied magnetic field, B z (0) is exactly the external applied field strength and the 
non- vanishing B z {k) comes solely from the non-linear current-phase response of the 
junction. 

Upon an integration over y, with the boundary condition A x (x, y — 0) = 
imposed by the gauge choice (|B.7[) . we obtain 



A x {x,y) =sgn(y) 



*> dfc 



(0)(C 



-HAl 



1 



2 B z {k) 



cos(kx + 6k) 



L-<*k\y\ _ A 



(B.13) 



/0+ 27T ' ■ ' Oik 

By using Eq. (|B.9|) and the boundary condition Jf(y —> ±00) = 0, we find that the 
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superconductor phase is given by 

2n r 
4>{x,y) = -sgn(y)— - A Lj B 2 (0)x 
*o L 

r°° Ak 2\B z {k)\sm{kx + 6 k ) 
From ip(x) = 2<f>(x, + ), we thus have that the phase difference across the junction is 



Aj^e-"^! 



1 



(B.14) 



tp(x) = - 



4ttA| 



xB z {0) 
Al 



dfc 2\B z (k)\a k 
2n k 



sin(fca: + 9k) 



(B.15) 



This result immediately provides a way to connect the phase difference and the 
magnetic-field distribution at the junction. 

Two important consequences can be drawn now. First, we make the observation: 

Ho 2en s h d 2 ip(x) f°° dfc 
~2 ~m : 



dx 2 



()■!■ 



2tt 



2 \B Z (k)\ hat sin (kx 



(B.16) 



In the long wavelength limit, fc <C 1/Al, the expansion a k = \/k 2 + 1/A 
±-(1 + ^X1 + 



.) leads to the relation 
dB z (x) efi n s h\ L d 2 ip(x) 



O 



d 4 (p(x) 



(B.17) 
(B.18) 

2 7 r W) J c (2A L ) aild 5 = 

, and using dimcnsionlcss variables ^ = x/Xj and r = ct/Xj, we arrive at 



9x to* <9x 2 \ 9x 4 

By inserting this relation into Eq. (|B.3|) , we obtain the equation of motion 

en n s hX L d 2 ip(x,t) . 

— ; ^2— = ^o>/c sin (^(x, 0) + ■ ' ' 

Then by defining the corresponding length scales Aj = 



jgo 



e/io(2A L +«) 

the dimcnsionless sine-Gordon equation 



d 2 d 2 d 

d^~d^~ a d^^ ^ (C ' T) = Sin ^ (C ' r) 



(B.19) 



Here, a and 7 are defined in the main text, and the Josephson penetration depth Aj 
recovers the given in the main text provided the effect of the insulator thickness I is 
taken into account. 

Second, the boundary condition on the phase difference can be inferred from 
Eq. (JB.15J) . The phase consists of two contributions: a term linearly increasing with 
x and an oscillating one. For a wide junction, i.e., w 3> Aj, the contribution of the 
oscillating term is negligible (on average) by comparison to the linear term. This 
translates into the condition 

V ( W ) - p(0) = -2tt%£ (B.20) 

where <I> c ff = (2X^)wB z (0) is the effective total flux threading through the junction 
area. Because B z (0) is the average applied magnetic field, the boundary condition for 
Eq. (|B.18J) is controlable externally. Again, the effective total flux yields an additional 
correction with 2Al — > 2Al + 1 once the effect of the insulator thickness I is taken into 
account. 

In the limit a,7 C 1, the phase difference profiles can be approximated by first 
neglecting those terms. The solution of sine-Gordon equation can be cast in term of 
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elliptic integrals (37J with the boundary condition (|B.20[) . In Fig. IB1[ we plot a typical 
phase profile for an effective flux density of $o/10Aj, as a function of the dimcnsionless 
length (. We observe that solitons (Josephson vortices) form, separated by lOAj. 

Appendix C. Fermionic states bound to a Josephson vortex in a 
topological superconductor 

Within the topological superconductor region, the fermionic states bound to a 
Josephson vortex can be understood via a simple edge-state coupling model |21j . 
The low-energy fermionic degrees of freedom along a wide Josephson junction can be 
described by the Hamiltonian 

H — w^ I dx (ipRd x ipii - tPlOx^l) +2im / dxcos(cp/2)ip R x/j L . (C.l) 

where ipu(L) lS the right(lcft) moving Majorana fermion operator, v^ is the velocity of 
the Majorana edge states and m is the tunnelling amplitude. The first term is the free 
Hamiltonian of a pair of counter-propagating Majorana edge states, and the second 
term accounts for the tunnelling amplitude, which depends on the superconducting 
phase difference ip across the junction. 

The Hamiltonian (|C.1|) satsifies the quantum mechanical supersymmetry that 
guarantees the presence of a Majorana zero-energy mode [21]. This zero- mode is 
responsible for the intrinsic entropy, (&b/2) In 2. To obtain the low-energy excitations, 
we first note that the the Josephson vortex profile, ip = 4 tan -1 exp(x/Z/), leads to 
cos(</?/2) = — ta,nh(x/lf). By linearizing the tunnelling term with the profile of tp, one 
can solve for the energy spectrum, to obtain 



s j = ± l^m n n = 012) ^ (C2) 

with a degeneracy g n = 1 for each state |43| . Hence, the minigap can be estimated as 



j _ 2fe^m / 2fe A T i , s. 

where we have approximate the tunnelling amplitude by m ~ Ati, appropriate for a 
transparent junction. From Sec. |Appendix A| we have that the entropy contributions 
due to the presence of the minigap states can be estimated from this minigap. 
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